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1.0  EXECUTIVE  SUMMARY 


1.1  Overview  of  the  Project 

Future  Air  Force  networks  will  face  the  challenge  of  providing  robust  data  and  circuit  services 
to  tens  or  hundreds  of  fixed  and  mobile  users  with  different  service  levels.  Some  of  the  service 
challenges  include  guaranteed  rates,  communication  over  difficult  channels,  hard  time-deadlines, 
reliable  message  delivery  over  unreliable  networks,  security,  and  policy-driven  resource  allocation. 

In  light  of  the  propagation  characteristics  as  well  as  the  new  challenges  of  the  aerial  layer 
communications,  the  project  “Cooperative  Routing  for  Dynamic  Aerial  Layer  Networks”  is  funded 
by  Air  Force  Research  Laboratory  (AFRL)  under  grant  number  FA8750-14-1-0077  to  develop  a 
novel  and  fundamental  framework  on  which  new  enabling  aerial  layer  communication 
technologies  and  protocols  can  be  designed  and  analyzed.  To  be  specific,  the  project  1)  focuses  on 
dynamic  spectrum  access  (DSA)  networks  which  efficiently  allow  different  wireless  nodes  to  co¬ 
exist  in  the  same  radio  spectrum;  2)  introduces  a  new  cooperative  aerial  layer  routing  strategy 
based  on  mutual-information  accumulation  (e.g.,  rateless  codes,  fountain  codes,  and  hybrid  ARQ); 
3)  develops  dynamic  network  resource  management  and  collaborative  routing  strategies  in  support 
of  aerial  layer  networking;  and  4)  demonstrates  the  effectiveness  of  the  strategies  introduced  using 
Ettus  universal  software  radio  peripheral  (USRP)  N210s.  Our  analytical  results  show  that  both  the 
centralized  and  distributed  algorithms  can  reduce  up  to  77%  of  the  end-to-end  delay  compared  to 
the  traditional  routing  strategies  in  DSA  networks;  demo  results  using  USRP  N210s  confirm 
anal5^ical  findings. 


1.2  List  of  People  Involved 

There  is  one  faculty  member  and  three  Ph.D.  students  involved  in  the  project,  which  are 
listed  in  the  following: 

■  Lingjia  Liu;  University  of  Kansas  (KU);  PI 

■  Hao  Chen;  University  of  Kansas  (KU);  Ph.D.  student  investigator 

■  Rachad  Atat;  University  of  Kansas  (KU):  Ph.D.  student  investigator 

■  Somayeh  (Susanna)  Mosleh;  University  of  Kansas  (KU):  Ph.D.  student  investigator 

It  is  important  to  note  that  both  Hao  and  Rachad  defended  their  Ph.D.  in  June  2017  with  help  from 
the  grant. 
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2.0  INTRODUCTION 


In  the  Air  Force’s  network-enabled  paradigm,  many  capabilities  will  be  generated  through,  and 
dependent  on,  the  integrated  efforts  of  multiple  components.  In  general,  robust  dynamic  tactical 
networks  can  support  10,000  to  500,000  communicating  devices  in  theater,  going  beyond 
traditional  tactical  networks,  oriented  to  the  voice  and  human  users,  and  beyond  traditional  ad  hoc 
networks.  Furthermore,  the  Air  Force  network  must  provide  robust  data  and  circuit  services  to  tens 
or  hundreds  of  fixed  and  mobile  users  with  different  service  levels.  Some  of  the  service  challenges 
include  guaranteed  rates,  communication  over  difficult  channels,  hard  time-deadlines,  reliable 
message  delivery  over  unreliable  networks,  security,  and  policy-driven  resource  allocation.  These 
characteristics  are  primarily  driven  by  defense  needs. 

In  response  to  challenges  associated  with  the  current  and  anticipated  Air  Force  network,  the 
concept  of  aerial  layer  networking  is  introduced.  The  aerial  layer  provides  a  new  horizon  for  aerial 
wireless  communications.  By  utilizing  aerial  layer  assets  as  communications  nodes  and  relays,  the 
Air  Force  and  its  partners  can  readily  amplify  and  augment  available  communications  assets  and 
enhance  warfighter  collaboration.  Using  the  aerial  layer,  communications  assets  are  not 
constrained  by  terrestrial  topology  and  the  limitations  hindering  ground  and  space-based 
communications.  Instead,  the  aerial  layer  enables  dynamic  network  topology  changes,  increases 
bandwidth  on-the-fly,  and,  additionally,  alternates  routing  for  existing  communications  networks 
to  reduce  enemy  threats. 

Dynamic  spectrum  access  (DSA),  also  known  as  Cognitive  Radio  Network  (CRN),  has  been 
regarded  as  one  of  the  key  techniques  to  improve  the  spectral  efficiency  of  a  wireless  network.  The 
basic  idea  of  DSA  is  to  allow  secondary  users  (SUs)  to  use  primary  user  (PU)  spectrum  under  the 
condition  that  SUs  will  not  cause  unacceptable  interference  to  PUs.  Several  approaches  have  been 
introduced  on  how  SUs  should  access  the  licensed  spectrum  such  as  spectrum  overlay 
[Liang2008Sensing]  [Rem2014Cad],  spectrum  underlay  [Sharma2008Channel]  or  cooperation 
[Su2010Active].  Overlay  approaches  enable  SUs  transmit  only  when  PUs  are  sensed  as  idle. 
Underlay  approaches  allow  SUs  to  coexist  with  PUs  under  the  constraint  that  SU’s  interference  to 
PUs  is  under  a  certain  threshold.  For  cooperation  approaches,  SUs  act  as  a  relay  to  collaborate 
with  the  PUs. 

In  underlay  DSA  networks,  SUs  and  PUs  are  transmitting  simultaneously  on  the  same  radio 
spectrum.  In  order  to  make  sure  that  the  interference  from  SUs  to  PUs  is  tolerable,  SUs  are  usually 
required  to  transmit  at  a  low  power  level  resulting  in  a  rather  limited  link  coverage.  To  address  this 
issue,  in  this  project,  we  introduce  cooperative  routing  strategies  among  SUs  to  extend  SU  network 
coverage.  To  be  specific,  we  consider  a  SU  network  that  employs  advanced  channel  coding 

2 


Approved  for  Public  Release;  Distribution  Unlimited. 


strategies  such  as  mutual-information  accumulation  at  the  physical  layer,  and  study  the  cooperative 
routing  and  resource  allocation  problems  associated  with  such  SU  networks. 

Rateless  codes  are  an  advanced  technique  compared  with  fixed-rate  codes  and  a  realistic 
method  to  realize  mutual-information  accumulation  in  the  physical  layer.  With  a  fixed-rate  code, 
the  transmitter  needs  to  design  codes  based  on  the  dynamic  channel  state  information  (CSI),  often 
incurring  resource  waste  or  too  many  errors.  To  guarantee  reliable  transmission,  the  data  length 
and  number  of  redundancy  bits  are  always  long,  decreasing  the  transmission  efficiency.  In  contrast, 
rateless  codes  allow  the  transmitter  to  generate  an  unlimited  number  of  encoded  packets,  such  that 
the  receiver  can  decode  the  data  after  receiving  a  sufficient  number  of  encoded  packets, 
irrespective  of  which  ones  it  has  received  [Shokrollahi2006Raptor].  In  addition,  rateless  codes  do 
not  require  knowledge  of  the  CSI  at  the  transmitter,  and  offer  greater  robustness,  reliability  and 
efficiency  as  compared  to  fixed-rate  codes  [Luby2002LT,  Erez2012Rateless]. 

Robustness  is  a  very  important  issue  in  a  communication  network,  especially  for  DSA 
networks.  The  network  should  continue  to  function  even  when  some  of  the  communication  nodes 
are  destroyed  or  turned  off.  This  means  that  nodes  in  the  network  could  be  able  to  work  in  the  peer- 
to-peer  (P2P)  or  device-to-device  (D2D)  mode  and  each  node  in  the  network  could  collaborate 
with  other  nodes  in  a  distributed  fashion.  Therefore,  after  characterizing  the  optimal  resource 
allocation  and  optimal  cooperative  route  for  one  source-destination  pair,  we  will  investigate 
distributed  routing  and  network  scheduling  algorithms,  and  extend  it  to  multiple  source-destination 
pairs.  The  introduced  algorithm  will  be  distributed  in  the  sense  that  there  is  no  centralized 
controller  that  has  to  have  all  information  about  the  whole  network.  When  designing  distributed 
algorithms,  the  objective  in  this  project  becomes  designing  the  distributed  algorithm  in  a 
systematic  way.  In  general,  distributed  algorithms  can  be  identified  through  designing  a  routing 
metric  by  experience.  Even  though  a  good  routing  metric  yields  a  reasonable  distributed  algorithm 
in  certain  scenarios,  it  is  far  from  systematic.  Throughout  this  project,  our  attention  focuses  on 
designing  the  distributed  algorithms  systematically,  that  is,  designing  the  distributed  algorithm 
based  on  the  decomposition  of  the  optimal  solution  from  the  centralized  problem. 

It  is  important  to  note  that  five  USRP  N2I0s  were  borrowed  from  AERL  to  perform  the 
hardware  demonstration  at  the  end  of  this  project.  Two  demonstrations  were  done  at  Technical 
Interchange  Meetings  (TIM)  in  June  2015  and  May  2016.  More  specifically,  cooperative  routing 
using  rateless  codes  (Raptor  Codes)  in  a  five-node  DSA  network  was  demonstrated  and  the 
performance  gain  of  the  hardware  demonstration  was  shown  to  be  in  strong  agreement  with 
theoretical  results. 
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The  project  was  organized  around  four  interconnected  research  thrusts: 

1)  Thrust  1:  spectrum  sensing  for  dynamic  spectrum  access  networks; 

2)  Thrust  2:  cooperative  routing  and  resource  allocation  based  on  rateless  coding  for  dynamic 
spectrum  access  networks; 

3)  Thrust  3:  distributed  network  resource  management  for  rateless  coding-based  cooperative 
dynamic  spectrum  access  networks; 

4)  Thrust  4:  USRP  demo  using  different  waveforms,  modulation  schemes,  and  multi-access 
strategies. 

In  the  following  sections,  we  describe  the  Methods,  Assumptions,  Procedures,  Results  and 
Discussions  for  each  of  the  aforementioned  contributions. 
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3.0  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 


The  project  can  be  summarized  into  four  interconnected  research  thrusts:  Thrust  1:  spectrum 
sensing  for  dynamic  spectrum  access  networks;  Thrust  2:  cooperative  routing  and  resource 
allocation  based  on  rateless  coding  for  dynamic  spectrum  access  networks;  Thrust  3:  distributed 
network  resource  management  for  rateless  coding-based  cooperative  dynamic  spectrum  access 
networks;  and  Thrust  4:  USRP  demo  using  different  waveforms,  modulation  schemes,  and  multi¬ 
access  strategies. 

The  following  technologies  are  adopted  in  the  project  to  advance  the  objective  in  each  thrust. 


Mutual  Information  Accumulation/Combination  Versus  Ener2v  Accumulation 

In  the  traditional  automatic  repeat  request  (ARQ)  transmission,  if  a  data  packet  is  corrupted 
and  cannot  be  decoded  error  free  at  the  receiver,  the  same  packet  will  be  requested  to  be 
retransmitted.  All  the  previous  receptions  of  the  packet  will  be  discarded  and  this  procedure  will 
repeat  until  the  packet  is  received  error  free.  Since  each  reception  of  the  packet  contains 
information  about  the  original  packet  even  when  it  is  not  decoded  error  free,  it  is  much  more 
efficient  for  the  receiver  to  keep  the  received  copy  instead  of  throwing  it  away.  This  is  the  main 
motivation  of  hybrid  ARQ.  There  are  two  classes  of  hybrid  ARQ:  chase  combining  and  incremental 
redundancy. 

•  Chase  combinins:  every  retransmission  contains  the  same  information.  The  receiver 
uses  maximum-ratio  combining  to  combine  the  received  data.  Because  all  transmissions  are 
identical,  chase  combining  can  be  seen  as  additional  repetition  coding.  One  could  think  of 
every  retransmission  as  adding  extra  energy  to  the  received  transmission.  Therefore,  this  mode 
can  also  be  regarded  as  energy  accumulation. 

•  Incremental  redundancy:  every  retransmission  contains  different  information  than  the 
previous  one.  Multiple  sets  of  coded  bits  are  generated,  each  representing  the  same  set  of 
information  bits.  The  retransmission  typically  uses  a  different  set  of  coded  bits  than  the 
previous  transmission,  with  different  redundancy  versions  generated  by  puncturing  the  output 
of  the  encoder.  Therefore,  at  every  retransmission  the  receiver  gains  extra  information.  This 
mode  can  also  be  regarded  as  mutual  information  accumulation/combination. 
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The  difference  between  energy  accumulation  and  mutual  information  accumulation/ 
combination  is  most  easily  understood  by  considering  binary  signaling  over  a  pair  of  independent 
erasure  channels.  Two  cooperating  transmitters  wish  to  transmit  a  common  message  to  a  single 
destination.  If  the  erasure  probabilities  are  both  pg,  and  both  transmitters  use  the  same  code,  then 
each  symbol  will  be  erased  with  probability  p| .  Therefore,  1  —  p|  novel  parity  symbols  are 
received,  on  average,  per  transmission.  If,  instead,  the  two  transmitters  use  different  codes,  on 
average  2(1  —  pg)  novel  parity  symbols  (which  exceeds  1  —  p|)  are  received  per  transmission. 
The  latter  is  mutual-information  accumulation/combination,  while  the  former  is  an  example  of 
energy  accumulation.  From  this  simple  example,  it  is  obvious  that  mutual  information 
accumulation/combination  is  strictly  better  than  energy  accumulation. 

Network  Resource  Allocation  under  Mutual  Information  Accumulation/Combination 

When  the  nodes  in  a  network  are  equipped  with  mutual  information  accumulation/combination, 
the  cooperative  communication  problem  becomes  a  “roM/mg”  problem,  i.e.,  identifying  which  of 
the  available  nodes  should  participate  in  the  transmission  and  what  system  resources  (time,  energy, 
and  spectral  bandwidth)  should  be  allocated  to  each.  It  is  clear  that  the  physical/medium-access- 
control  layer  technique  employed  strongly  influences  the  optimum  route. 

There  has  been  little  prior  work  investigating  routing  in  networks  consisting  of  nodes  using 
mutual  information  accumulation/combination.  Mutual  information  accumulation/combination  is 
investigated  in  [Molisch2007Performance],  but  the  analysis  therein  assumes  network  “flooding”, 
i.e.,  all  nodes  transmit  all  the  time;  this  is  not  an  optimum  use  of  energy.  In  [Zhao2005 Practical]  a 
heuristic  algorithm  for  relaying  information  with  chase  combining  hybrid  ARQ  over  time  is 
derived.  In  contrast  to  the  technique  we  are  going  investigate  in  the  proposal,  [Zhao2005Practical] 
assumes  that  when  relay  nodes  transmit  simultaneously,  they  send  out  the  same  signal. 

When  considering  the  general  task  of  optimizing  route  and  resource  allocations  for  a 
communication  network  with  arbitrary  attenuation  between  nodes,  our  strategy  is  to  introduce  first 
the  idea  of  “transmission  order”  as  illustrated  in  [Draper20 11  Cooperative].  In  general, 
“transmission  order”  is  the  order  in  which  the  nodes  are  allowed  to  come  on-line  as  transmitters. 
We  can  think  of  the  transmission  order  as  the  route  used  by  the  cooperative  scheme.  Since  a  node 
cannot  transmit  until  it  has  decoded  the  message,  a  node’s  position  in  the  transmission  order  puts 
constraints  on  the  resources  allocated  to  transmitters  prior  to  it  in  the  order.  We  then  iterate  between 
two  sub-problems: 

1)  First,  for  the  given  transmission  order,  we  determine  the  optimum  transmission  parameters. 

This  resource  allocation  problem  turns  out  to  be  a  linear  program  (LP). 

2)  Second,  based  on  the  solution  of  the  LP  we  revise  the  transmission  order. 
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Distributed  Network  Resource  Allocation  for  Robustness 

Robustness  is  a  very  important  issue  in  a  communication  network,  especially  for  aerial  layer 
networks.  The  network  should  continue  to  function  even  when  some  of  the  communication  nodes 
are  destroyed  or  turned  off.  This  means  that  nodes  in  the  network  could  be  able  to  work  in  the  P2P 
or  D2D  mode  and  each  node  in  the  network  could  collaborate  with  other  nodes  in  a  distributed 
fashion.  Therefore,  after  characterizing  the  optimal  resource  allocation  and  optimal  cooperative 
route  for  one  source-destination  pair,  we  will  investigate  distributed  routing  and  network 
scheduling  algorithms,  and  extend  it  to  multiple  source-destination  pairs.  The  introduced  algorithm 
will  be  distributed  in  the  sense  that  there  is  no  centralized  controller  that  has  to  have  all  information 
about  the  whole  network. 

When  designing  distributed  algorithms,  the  objective  in  this  proposal  is  to  design  the 
distributed  algorithm  in  a  systematic  way.  In  general,  distributed  algorithms  can  be  identified 
through  designing  a  routing  metric  by  experience.  Even  though  a  good  routing  metric  yields  a 
reasonable  distributed  algorithm  in  certain  scenarios,  it  is  far  from  systematic.  Throughout  this 
project,  our  attention  focuses  on  designing  the  distributed  algorithms  systematically,  that  is, 
designing  the  distributed  algorithm  based  on  the  decomposition  of  the  optimal  solution  from  the 
centralized  problem. 


3.1  Spectrum  Sensing  for  Dynamic  Spectrum  Access  Networks 

A  sensing-based  spectrum  sharing  (SBSS)  system  [Kang2009Sensing]  enables  SUs  to  adapt 
their  transmit  power  according  to  their  detection  of  the  PUs’  status  as  well  as  the  channel  conditions. 
Compared  with  the  traditional  spectrum  sharing,  SBSS  is  more  flexible  and  dynamic,  allocate 
system  resources  to  maximize  overall  system  throughput. 

3.1.1  The  protection  of  PUs  and  the  sensing  throughput  trade-off  in  SBSS 

There  are  two  key  issues  in  SBSS  system:  the  protection  of  PUs  and  the  sensing-throughput 
trade-off.  To  protect  the  PU,  an  interference  power  constraint  or  interference  temperature 
constraint  is  imposed  on  SUs  in  spectrum  underlay,  while  in  spectrum  overlay,  a  high  detection 
probability  is  required  in  order  to  maintain  PUs’  spectrum  priority.  Therefore,  interference  power 
constraints  are  considered  in  underlay  systems  and  miss  detection  constraints  are  considered  in 
overlay  systems.  Neither  of  these  constraints  can  directly  measure  the  Quality-of-Service  (QoS)  at 
which  cellular  links’  performance  is  affected  by  D2D  users.  As  a  result,  a  rate  loss  constraint  (RLC) 
is  introduced  in  [Kang20100ptimal]  as  a  new  metric  which  can  be  directly  related  to  PUs’  QoS. 
This  new  criterion  is  shown  to  provide  improved  capacity  performance  over  conventional  criteria. 
We  adopt  a  rate  loss  constraint  in  our  SBSS  system. 
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Another  issue  in  SBSS  systems  is  the  sensing-throughput  trade-off  which  is  the  balance 
between  the  protection  of  PU’s  performance  and  the  maximization  of  SUs’  throughput.  The 
sensing  module  should  select  the  best  spectrum  resource  to  sense  and  optimize  the  sensing  time  in 
order  to  maximize  SUs’  throughput. 

In  [Liang2008Sensing]  and  [Pei2009How],  the  sensing  throughput  trade-off  in  spectrum 
overlay  under  the  detection  probability  constraint  is  investigated.  In  [Hamdi2009Power],  an 
algorithm  to  find  the  optimal  sensing  time  and  power  allocation  for  spectrum  overlay  using 
interference  power  constraint  is  introduced.  In  [Stotas20 11  Optimal],  the  problem  of  optimizing 
sensing  time  and  power  allocation  in  SBSS  cognitive  radio  networks  under  an  interference  power 
constraint  is  studied.  In  [Shi2012Joint],  an  optimal  power  and  sub-carrier  allocation  strategy  to 
maximize  SUs’  throughput  subject  to  SUs’  QoS  constraint  as  well  as  PU’s  RLC  constraint  is 
introduced  for  multi-carrier  systems. 

3.1.2  System  model  for  SBSS 

The  SBSS  system  consists  of  four  communication  nodes:  primary  user  transmitter  (PUT), 
primary  user  receiver  (PUR),  secondary  user  transmitter  (SUT),  and  secondary  user  receiver 
(SUR),  as  shown  in  Fig.  1.  As  in  most  wireless  systems,  a  frame-based  scheme  is  employed  in  the 
SBSS  system  where  the  frame  structure  consists  of  a  sensing  slot  x  and  a  data  transmission  slot  T 
-  X.  During  the  sensing  slot,  energy  detection  is  adopted  at  the  SUT  to  detect  whether  the  jth  sub¬ 
carrier  is  occupied  by  a  PU.  The  two  hypotheses  of  spectrum  sensing  can  be  summarized  as: 

=  zim)  (1) 

y(m)  =  m)  +  z(m)  (2) 

where  hp^^  is  the  channel  coefficient  between  PUT  and  SUT,  x(J,m)  is  the  mth  sample  of  the 
transmitted  signal  in  the  yth  sub-carrier,  and  z(m)  is  the  mth  sample  of  the  noise.  Denoting  the 
transmit  power  of  PU  as  Pp,  we  assume  that  both  the  received  signal  and  noise  are  zero-mean 

circularly  symmetric  complex  gaussian  with  variance  =  gpsfPp  and  =  rj,  respectively.  The 
test  statistic  for  an  energy  detector  is  given  by 


1 

r(y)  =  — ^  |y(m)p 


m=l 


(3) 


where  fs  is  the  sampling  frequency.  The  false  alarm  probability  and  detection  probability  in  the  ;th 
sub-carrier  can  be  written  as 
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(4) 


p/(£,t)  =  Pr(r(y)  >  e\Ho)  =  -  l) 

pj^s,  t)  =  Pr(r(y)  >  e\H,)  =  q[[^-  W  “  l)  j  (5) 

where  is  the  ratio  between  power  of  PU’s  signals  received  at  SU  and  noise,  Q{-^,  is  the 
complementary  distribution  function  of  the  standard  Gaussian.  To  simplify  the  analysis  in  the 
following,  the  sensing  threshold  is  fixed  as  e  =  (1  +  /2)(7^,  which  ensures  that  >  Pf. 


Figure  1:  The  SBSS  system  model. 

Let  P(Ho)  denote  the  probability  that  PU  is  idle  and  P(//i)  denote  the  probability  that  PU  is 
busy.  If  SUT  senses  PU  idle  in  sub-carrier  j,  it  will  transmit  at  a  power  level  of  P^-,  otherwise,  it 
will  transmit  at  a  power  level  of  P/.  Let  gpp,gps>gip>gis  represent  the  channel  gains  of  the  ;th 
sub-carrier  between  PUT  and  PUR,  PUT  and  SUR,  SUT  and  PUR,  SUT  and  SUR,  respectively. 
Based  on  PU’s  activity  and  SU’s  sensing  results,  the  throughput  of  SU  can  be  computed  according 
to  the  following  four  scenarios: 

1)  Scenario  1:  SU  successfully  detects  that  PU  is  idle.  The  probability  of  this  case  is  = 
P(Hq)(1  —  P^ (t)),  SU’s  achievable  rate  can  be  expressed  as  r^Q  =  log2(l  +  P^gis/g)- 

2)  Scenario  2:  SU  makes  a  false  alarm.  The  probability  of  this  case  is  a{  =  P(Hq)PI  (t).  SU’s 
achievable  throughput  can  be  expressed  as  =  log2(l  +  P/ gis/g) 

3)  Scenario  3:  SU  fails  to  detect  PU’s  presence.  The  probability  of  this  case  is  Pq  = 
P(//i)(l  —  P^ (t)).  SU’s  achievable  rate  can  be  expressed  as  t/q  =  log2(l  +  P^gis/ig  + 
Ppgis))- 

4)  Scenario  4:  SU  successfully  detects  PU’s  presence.  The  probability  of  this  case  is  Pi  = 

P(//i)P^(t).  SU’s  achievable  rate  can  be  expressed  as  +  P^ gl^/ir]  + 

P^gis)) 
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3.1.3  Power  Allocation  and  Sensing  Time  Optimization  Algorithm 


In  this  section,  we  study  the  optimal  sensing  time  and  power  allocation  strategy  for  SBSS 
networks  with  PU’s  RLC.  Specifically,  in  the  single-carrier  case,  we  will  characterize  the  optimal 
system  operation  and  compare  the  performance  (SU’s  optimal  throughput  and  optimal  sensing 
time)  of  SBSS  with  that  of  a  spectrum  overlay  network.  In  the  multi-carrier  case,  we  will  introduce 
an  efficient  algorithm  to  find  the  optimal  sensing  time  and  power  allocation. 

As  the  first  step,  we  focus  on  the  single-carrier  model  where  the  channel  index  j  is  omitted.  In 
this  scenario,  the  SU’s  average  throughput  during  each  frame  can  be  written  as: 

T-t 

fir,  Pq,  Pi)  =  Kroo  +  -I-  (6) 


where  Tqq,  r^o,  Ui  SU’s  achievable  rate  when  SU  successfully  detects  that  PU  is  idle,  SU 
makes  a  false  alarm,  SU  fails  to  detect  PU’s  presence,  and  SU  successfully  detects  PU’s  presence, 
respectively. 

Denoting  6  as  the  percentage  of  rate  loss  that  PU  can  tolerate  and  C  =  T  log2  (1  +  igppPp) /??) 
as  PU’s  data  rate  without  SU’s  interference,  to  protect  PU’s  performance,  SU  should  guarantee 
that  PU’s  average  (over  sensing  results)  achievable  rate  should  not  be  less  than  (1  —  6)C,  which 
can  be  expressed  in  the  following  equation: 


hl(T,Po,Pi)  =  (r-T) 


(1  -  Pa)  log2 


9  sp  Pq  ) 

V  d"  9spPo) 


+  Pd  l0g2 


+  ^C-(1-5)C>0 


9 sp  Pq  \ 
V  +  9spPi) 


(7) 


In  order  to  meet  the  power  budget  requirement  at  the  SU,  an  average  power  constraint  (over 
all  sensing  results)  should  be  considered,  which  can  be  written  as: 

T -T 

h2(j,Po’Pi)  =  Ptot - j —  [(^0  +  Po)Po  +  (^1  +  ^  0  (8) 


where  Ptot  is  SU’s  average  power  budget. 

Accordingly,  in  the  single-carrier  SBSS  system,  to  maximize  SU’s  throughput  under  PU’s 
rate  loss  constraint  and  SU’s  power  constraint,  an  optimization  problem  OPl  can  be  formulated  in 
the  following: 


max  /(t,Po,Pi)  (9) 

{t,Po,Pi} 

s.t.:  /Ii(t,Po,Pi)  >  0,h2(T,Po,Pi)  >  O.Pq  >  0,P^  >  0,T  >  T  >  0 
In  the  multi-carrier  case,  let  J  denote  the  total  number  of  sub-carriers  in  the  system. 
Furthermore,  we  assume  that  both  PU  and  SU  can  occupy  several  channels  concurrently  based  on 
the  channel  conditions  and  sensing  results.  Accordingly,  SU’s  throughput  optimization  problem  in 
multi-carrier  case  under  PU’s  rate  loss  constraint  and  SU’s  transmission  power  constraint  can  be 
formulated  as: 
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(10) 


T  ^ 

^  +  Po^Io  +  «in/i  +  /^i^A] 

’  ^  j=i 

s.t:g^iT,Po,Pi)  >  0,g2i'r,Po,Pi)  >  0 
where  Pq  =  (^Pq.Pq,  ...,Pq),Pi  =  ...,P^).  5'i(t,  Pq,  Pi)  represents  PU’s  RLC  and 

Qiij’  Po>  Pi)  represents  SU’s  average  (over  sensing  results)  power  constraint. 

The  solution  to  find  the  optimal  sensing  time  and  power  allocation  in  single-carrier  case  is 
summarized  in  the  following  algorithm: 

Algorithm  1:  Find  optimal  sensing  time  and  power 
allocation  in  single-carrier  SBSS 

Input:  System  parameters:  T,  S,  P{Ho),  Ptot^  /«> 

9spi  Ops')  Opp-!  Pp} '■> 

Output:  Optimal  sensing  time  r  and  Power  allocation 
Po,Pi\ 

1  for  r  ^  0  to  r  do 

2  calculate  Si  by  (10); 

3  calculate  S2  by  (13); 

4  calculate  Ss  by  (14); 

5  find  the  optimal  power  under  r,  (Po(r),  Pi(t)),  by 

_  (15); 

6  T*  =  argmaxT>r>o/('r,Po(T),Pi(T)); 

7  {P^,P,*)  =  {Po{T*),Pi{T*)y, 

8  return  t*,{P^,P^). 

The  solution  to  find  the  optimal  sensing  time  and  power  allocation  in  multi-carrier  case  is 
summarized  in  the  following  algorithm: 
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Algorithm  2:  Find  optimal  sensing  time  and  power 
allocation  in  multi-carrier  SBSS 


Input:  System  parameters:  T,  5,  P{Ho),  r),  Ptot,  fs^ 
{9L^9ip,9i.„9^p,Pi,j  = 

Output:  Optimal  sensing  time  r*  and  Power  allocation 
Po*,Pi* 

1  Initialization: 

Po(0),Pi(0),0  <  1/  <  1,0  G  (0,  l),r/ 

c(0)>  0,  /?(0)  -(-  //  -I-  1; 

2  for  T  <—  0  to  T  do 

3  while  P{k  —  1)  >  iz  do 

4  Start  from  (Po(fc  —  l),Pi(A;  —  1)),  solve  the 
unconstrained  optimization  problem  (18)  using 
BFGS  method  and  get  the  minimal  point  as 
(Po(fc),Pi(fc)); 

5  if  0{k)  >  6/3{k  —  1)  then 

6  |_  c(/c  -I- 1)  r]c{ky, 

7  else 

8  L  ^ 

9  update  1.1  j  by  (19); 

10  k  ^ —  A:  -|-  1; 

11  L  (Po*(T),Pi*(r))  ^  (Po(A-),Pi(fc)); 

12  r*  ■«- argmax7’>T>o/(T,Po*(T),Pi*(r)); 

13  (Po*,Pi*)^(Po*(r*),Pi*(r*)); 

14  return  t*,Po*,Pi*; 
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3.2  Rateless  Coding  based  Cooperative  Dynamic  Spectrum  Access  Networks 

As  introduced  in  the  former  sections,  mutual-information  accumulation  is  an  enhanced 
technique  to  increase  the  performance  of  the  dynamic  spectrum  access  network  which  can  be 
realized  through  rateless  codes.  Rateless  codes  allow  the  transmitter  to  generate  an  unlimited 
number  of  encoded  packets,  such  that  the  receiver  can  decode  the  data  after  receiving  a  sufficient 
number  of  encoded  packets,  irrespective  of  which  ones  it  has  received  [Shokrollahi2006Raptor]. 
In  addition,  rateless  codes  do  not  require  knowledge  of  the  channel  state  information  (CSI)  at  the 
transmitter,  and  offer  robustness,  reliability  and  efficiency  as  compared  to  fixed-rate  codes 
[Luby2002LT,  Erez2012Rateless]. 

3.2.1  Improving  the  throughput  of  CRN/DSA  with  mutual-information  accumulation 

In  DSA  networks,  unlicensed  users  or  secondary  users  (SUs)  periodically  sense  the  spectrum 
to  detect  any  licensed  or  primary  user  (PU)  activity  before  transmission.  A  SU  needs  to  maximize 
its  opportunistic  transmissions  while  being  “transparent”  to  PU’s  activities.  In  practice,  a  primary 
receiver  (PR)  may  fail  to  decode  the  packets  from  a  primary  transmitter  (PT)  due  to  the  interference 
of  SUs  or  “too  weak”  channel  condition.  To  increase  SU  transmission  opportunities  and  increase 
PU  throughput,  the  SU  can  relay  PU  packets  that  are  unsuccessfully  received  at  the  PR.  This  would 
help  free  up  PU’s  queue  faster  so  SU  can  have  more  channel  access  [Simeone2007Stable, 
Zhang2013On,  Chen2015Optimal].  To  enhance  the  performance  of  PUs,  mutual-information 
accumulation  technique  can  be  used  at  the  PR. 

Primary  packet 


Figure  2:  A  general  system  model  that  a  pair  of  SUs  co-exit  with  a  pair  of  PUs. 

Fig.  2  shows  a  CRN  where  a  pair  of  SUs  coexit  with  a  pair  of  PUs  [Devroye2006 Achievable]. 
ST  senses  PT’s  activity  at  the  beginning  of  the  time  slot,  if  detected  idle,  ST  will  transmit  a  packet, 
provided  it  has  at  least  one  in  its  queue.  PT  and  ST  are  equipped  with  queues  of  infinite  length, 
with  ST  having  two  queues:  one  holding  its  own  packets,  and  the  other  holding  PT  packets  that 
were  unsuccessfully  received  by  PR. 

Independent  and  stationary  traffic  arrival  processes  are  considered  for  primary  and  secondary 
queues,  with  arrival  rates  Aj,  and  departure  rates  jUj  packets/sec,  where  i  reads  “P”  for  Qp  (holding 
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PT  packets),  “S”  for  Qs  (holding  ST  own  packets),  and  “PS”,  “SP”  for  Qps  (holding  PT  packets  to 
be  relayed  by  ST).  PT  transmits  with  normalized  power  Pp  =  1,  while  ST  selects  its  transmission 
power  P5  <  1.  A  packet  is  successfully  received  at  the  intended  destination  if  the  SNR  is  above  a 
threshold  13 1,  where  i  reads  “P”  for  primary,  and  “S”  for  secondary.  The  probability  of  an  outage 
event  is 


A_ 

rfi 


(11) 


Due  to  fading  and  path  loss,  ST  may  misdetect  PT’s  transmission  with  probability  p^,  causing 
undesirable  interference  to  the  primary  link.  Also  ST  might  detect  the  idle  slot  as  busy  with  a  false 
alarm  probability  Pf.  The  departure  rate  Xp(t)  of  the  primary  transmitter,  can  be  defined  as  the 
number  of  primary  packets  that  are  successfully  received  at  either  PR  or  ST  at  time  t.  Then 


E[X,(0]  =  E  {1  {o,(0  u  0;(/)})  = 


^  ^  ^  Pp  Pp  ^ 

e  +e  -e 


+  p 


/ 


rp 


(12) 


Op  (t)  and  Op  (t)  denote  the  events  that  PT  successfully  transmits  without  interference  from  ST 
and  under  ST  interference,  respectively. 


The  departure  rate  A5(t)  of  the  secondary  queue  (25  is  defined  as  the  number  of  secondary 
packets  that  are  successfully  transmitted  to  the  intended  SR  at  time  t,  and 


E[X,(()]  =  E{l{(4{()nO,{())u(.4j(()nOi(())}}  =  /i,(F,,e) 


A. 


e  ^"’'*(l-e)  +  p^ 


(l-€) 


PsYs^ 


zK 

s 


Ps7 s  PR 


(13) 


where  As(t)  and  denote  the  events  that  slot  t  is  available  for  transmission  by  ST  when 

correctly  sensed  as  idle  when  misdetected  as  idle,  respectively.  Osit)  and  denote  the  events 
of  successful  transmission  by  ST  under  no  interference  from  PT  and  under  PT  interference, 
respectively,  e  is  the  probability  that  ST  transmits  PT  packet  when  the  slot  is  sensed  idle 
[Elazzouni20 1 5Fullduplex] . 


Similarly,  X^p  (t)  is  defined  as  the  number  of  successful  reception  by  PR  from  ST  without 
interference  from  PT  and  under  PT  interference,  respectively.  Then 
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E[X,/0]  =  e{1{(^/0  n  O^/O)  u  (^'(0  n  o;(0)}}  = ,e) 

=  (\-Pf) 
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where  Osp  (t)  and  O^p  (t)  denote  the  events  of  successful  reception  by  PR  from  ST  with/without 
interference  from  PT,  respectively.  The  arrival  rate  of  Qps  is  defined  as  the  number  of  primary 
packets  that  are  to  be  relayed  by  ST,  then  the  average  arrival  rate  of  Qps  can  be  written  as 
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Xps  <  PspiPsi  should  be  guaranteed  to  keep  the  stability  of  the  queue  Qps(t),  then 
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The  above  optimal  e  to  maximize  the  throughput  of  the  CRN  is  valid  when  mutual-information 
accumulation  is  not  used.  When  mutual-information  accumulation  is  used,  everything  remains 
valid  except  the  probability  of  Osp{t).  In  mutual-information  accumulation,  nodes  accumulate 
information  as  bits/Hz,  where  Cij  is  the  channel  capacity  between  nodes  i  and  j,  then  the 
probability  of  0sp(t)  can  be  modified  as 
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where 
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Also,  to  guarantee  the  stability  of  the  queue  Ops(t),  e  >  ^ps/ksiPs)  should  be  met.  Note  that 
this  equation  can  be  solved  numerically. 


3.2.2  Cooperative  Routing  for  Underlay  CRNs  using  mutual-information  accumulation 

In  CRNs,  SUs  have  to  dynamically  control  their  transmit  powers  so  that  the  interference  to 
PUs  is  tolerable.  Under  this  constraint,  SUs’  data  link  usually  suffers  from  either  high  error  rate 
with  limited  transmission  range  or  long  end-to-end  delay  cause  by  multi-hop  transmissions.  To 
address  this  issue,  a  cooperative  routing  strategy  among  SUs  that  employs  advanced  channel 
coding  strategies  such  as  mutual-information  accumulation  is  applied  to  extend  SU  networks’ 
coverage  [Palanki2004Rateless,  Castura2007Rateless].  The  objective  of  the  SU  is  to  minimize  the 
end-to-end  delay  under  sum-bandwidth,  sum-energy  and  PU’s  interference  power  constraint. 


PUT  t  ■ 


PUR  r 


Figure  3:  System  model  of  Underlay  CRN  with  a  pair  of  PUs  and  N  SUs. 

Fig.  3  shows  a  underlay  CRN  containing  1  PU  transmitter,  1  PU  receiver,  1  SU  source,  1  SU 
destination  and  N  —  2  SU  relays.  The  spectrum  bandwidth  of  the  subband  is  denoted  as  Wp,  and 
PU’s  interference  power  threshold  is  Ip.  The  SU  network  needs  to  transmit  a  data  packet  with  an 
entropy  of  B  bits  from  the  SU  source  to  SU  destination  with  the  help  of  SU  relays.  Depending  on 
the  routing  and  resource  allocation  strategies,  all  SU  nodes  may  actively  help  forwarding  the 
packet  or  keep  silent.  Let  P5  be  the  transmit  power  spectral  density  of  SUs.  The  channel  power 
gain  between  node  i  and  j  is  hij .  Cij  is  the  achievable  transmission  rate  per  degree  of  freedom 
between  node  i  and  j 
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(19) 


where  Pp  is  the  transmit  power  spectral  density  of  PU,  f  is  the  SNR  gap  between  the  channel 
capacity  and  a  practical  modulation  and  coding  scheme. 


Let  L  be  the  total  number  of  nodes  in  the  transmission  order.  Aj=  Tj  —  Ti_i  is  the  inter-user 
delay  between  user  i  and  user  i-1.  Ti  is  the  time  at  which  user  i  decodes  the  message,  and  can 
be  regarded  as  the  length  of  the  ith  time  slot.  For  a  given  route  order,  the  objective  of  the  SU 
network  is  to  minimize  end-to-end  transmission  delay  Ti\ 


0P1:7;  =  2;A, 

1=1 


The  minimization  problem  is  constrained  by  the  following  constraints: 


k-l  k 

t=0  j=i+\ 


>B,  k  =  l,---,L. 


(20) 


(21) 


where  Aij  denotes  the  time-bandwidth  product  assigned  to  the  ith  user  in  the  jth  time  slot.  This 
equation  means  that  the  accumulated  mutual  information  should  exceed  the  entropy  after  L 
transmissions. 


1=0  y=i+i 


(22) 


this  constraint  means  the  total  energy  consumption  for  delivering  the  message  throughout  the  SU 
network  should  not  exceed  the  energy  budget  Ej  [Neely2005Dynamic]. 


7-1 


(23) 


this  constraint  suggests  that  all  the  nodes  in  the  SU  network  share  a  sum-bandwidth  of  Wp. 


1=0  (24) 

To  protect  the  PU’s  QoS,  SUs  should  make  sure  that  their  total  interference  powers  to  the 
PU  receiver  should  not  exceed  Ip.  hi^  is  the  channel  power  gain  from  the  /th  SU  to  the  PU  receiver. 
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In  fact,  OPl  is  NP-complete,  but  if  the  transmission  order  is  determined,  then  the  problem  is  easy 
to  solve.  Let  =  [Ai,  A2, ,  A^,  A*q^,  A*q2  - ,  ^i2>  - .  ^(l-i)l]  be  the  optimal  solution  of  OPl 

for  a  given  route/transmission  order.  Denote  the  corresponding  optimal  end-to-end  delay  as 


(25) 


Then  the  following  theorem  holds: 

•  Theorem  1:  if  A-  =  0,  let  Tl*  be  the  optimal  end-to-end  delay  of  the  “swapped”  route/ 
transmission  order 

[0, 2, /,/-l,/-l-l,---,Z]  if  /<  Z-1 
[O,-- -,1-2,1]  if i  =  L. 

then  Tl*  <  T^. 

•  Theorem  2:  If  user  i’s  transmission  time  is  0,  let  Tl*  denote  the  optimal  end-to-end  delay  of 
the  “deleted”  transmission  order: 

—  +  if  /<L-1 

thenrr  <  Tl. 

From  the  above  two  theorems,  we  know  that  for  any  route/transmission  order,  the  “sweeping” 
and  “deleting”  operation  ean  be  used  to  obtain  a  better  route/transmission  order. 

Algorithm  1  shows  the  detailed  iterative  route  optimization  and  resource  alloeation  algorithm. 
In  every  repeat,  the  OPl  problem  is  solved,  the  iteration  terminates  once  an  order  with  A^  >  0,  and 
tj  >  0  for  all  i  is  obtained. 
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Algorithm  1:  Centralized  cooperative  routing  and  resource 
allocation  algorithm 
Input:  System  parameters: 

Ps,  Pp,  a,  It^  Et^  |1  ^  hj  ^ 
iV},  {hir\l  <i<N-l},  {hpi\2  <  i  <  A^}; 
Output:  Route  from  the  SU  source  to  the  SU  destination 
and  the  corresponding  resource  allocation  Aij ; 

1  Initialization:  Apply  Dijkstra  algorithm  to  find  the  optimal 

path  from  the  source  to  the  destination  with  the  minimal 
interference  to  the  PU  receiver;  take  this  minimal 
interference  path  as  the  initial  transmission  order; 

2  repeat 

3  solve  the  linear  program  OPl  to  get  the  resource 
allocation  solution  for  current  transmission  order; 

4  (Based  on  Theorem  1  swap  the  transmission  order); 

5  if  =  0  and  A^_i  7^  0,  z  =  1, . . . ,  L  then 

6  |_  swap  the  positions  of  the  two  nodes  in  the  order; 

7  if  Al  =  0  then 

8  drop  node  L  —  1  from  the  order; 

9  (Based  on  Theorem  2  delete  the  silent  node); 

10  if  ti  =  0  then 

11  drop  node  i  from  the  order; 

12  until  An  order  with  A^  >  0  and  >  0  for  all  i  is 
obtained'. 


3.3  Distribute  Resource  Allocation 

Centrally  determining  resource  allocation  for  all  candidate  transmission  nodes  in  a  large-scale 
network,  if  even  possible,  is  an  exceedingly  expensive  operation.  A  potential  solution  is  to  perform 
distributed  resource  allocation  in  which  the  resource  allocation  is  determined  locally.  Compared 
with  the  centralized  resource  allocation,  the  distributed  scheme  is  a  quite  cost-efficient,  flexible, 
and  robust,  and  is  widely  used  in  many  research  fields  such  as  massive  machine  type 
communication  (MTC)  and  cooperative  routing  in  multi-hop  cognitive  radio  networks  (CRNs). 

3.3.1  Cooperative  Retransmission  for  massive  MTC 

In  massive  MTC  networks,  multiple  machine  type  devices  (MTDs)  need  to  transmit  packets  to 
their  base  stations  (BS)  [Andrews2014What].  The  wireless  connection  between  the  MTD  and  the 
BS  is  unreliable,  however,  due  to  the  interference  caused  by  uncoordinated  access  of  other  MTDs 
[Dhillon2015Wide].  Worse  yet  is  that  the  probability  that  the  retransmission  from  the  outage 
source  MTD  will  fail  again  due  to  the  spatiotemporally  correlated  interference  is  high 
[Ganti2012Spatial,  Nigam2015Spatiotemporal,  Haenggi2013Di versify].  It  is  impossible  for  the 
BS  to  centrally  perform  a  perfect  schedule  to  all  MTDs  to  avoid  such  interference.  A  potential 
solution  is  to  apply  a  distributed  location-based  cooperative  strategy  to  improve  the  performance 
of  massive  MTC  networks. 
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Figure  4:  System  model  of  the  distributed  cooperative  retransmission  strategy. 


Fig.  4  shows  the  system  model  of  the  cooperative  strategy.  In  the  cooperative  strategy,  an 
inactive  MTD  is  selected  as  a  relay  to  forward  the  packets  of  the  outage  source  MTD  in  the  first 
transmission.  The  relay  MTD  is  selected  based  on  the  following  two  requirements: 

1)  The  inactive  relay  MTD  has  successfully  decoded  the  packet  in  the  first  transmission; 

2)  The  relay  MTD  is  located  within  a  circular  area  Ay  around  the  BS  with  a  radius  of  Ry  to 
make  sure  the  relayed  packet  can  be  decoded  by  the  BS  with  a  high  probability. 

Assume  the  locations  MTDs  are  modeled  as  a  homogeneous  Poisson  point  process  (HPPP) 
4>  with  a  density  of  A  [Dhillon2014Fundamentals].  MTDs  transmit  collected  data  to  BS  with 
probability  p.  The  typical  BS  fails  to  decode  a  packet  if  the  received  SIR  is  less  than  a  threshold 
T,  and  all  MTDs  transmit  with  the  same  power  P^.  Therefore,  the  received  power  at  the  typical  BS 
from  a  MTD  at  a  distance  of  z  is  Py  =  where  Gs4_i~exp(l)  is  the  small  scale  fading 

gain  between  the  typical  MTD  s  and  the  typical  BS  d  during  the  first  transmission. 

Under  the  designed  cooperative  strategy,  the  typical  BS  is  outage  if:  1)  there  is  a  potential 
relay  selected  but  the  typical  BS  fails  to  decode  the  packet  after  the  retransmission  from  the 
potential  relay;  or  2)  there  is  no  potential  relay  and  the  typical  BS  fails  to  decode  the  packet  after 
the  retransmission  from  the  typical  MTD.  Thus,  the  outage  probability  of  the  designed  cooperative 
strategy  can  be  expressed  as 

pco  ^  pj'g)  -  P(A  n  B  n  F)  -  P(A^  n  B  n  F).  (26) 

In  the  equation,  A  is  the  event  that  there  is  at  least  one  potential  relay.  B  denotes  the  event  that  the 
first  transmission  from  the  typical  MTD  fails.  E  means  the  transmission  from  the  selected  relay 
succeeds,  and  F  indicates  that  the  retransmission  from  the  typical  MTD  succeeds.  A^  denotes  the 
complement  of  A  and  P(B)  =  1  —  e~‘^^  is  the  outage  probability  of  the  first  transmission.  After 
the  derivation,  the  formulation  of  B(A  n  B  n  E)  and  P{A^  n  B  n  F)  can  be  listed  as  follows 
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where 


^(3^) = (^2(7 + 17)  -  820, y)X 

t=i  >1  (29) 

Therefore,  in  a  massive  MTC  network,  the  outage  probability  for  a  typical  MTD  link  under  the 
distributed  cooperative  strategy  can  be  expressed  as 


oo  it+1 


^ ^ (  l)^g  <2  (A:  +  l) ^  pA(g(j>(i  p)x(i)  p(g(j+i)-g(j))) 
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(30) 


3.3.2  Distributed  Cooperative  Routing  and  Resource  Allocation  in  Underlay  CRNs 

In  Section  3.2.2,  the  cooperative  routing  for  underlay  CRNs  using  mutual-information 
accumulation  is  described,  but  the  optimal  cooperative  routing  and  resource  allocation  algorithm 
is  proposed  centrally.  The  operation  is  inefficient  in  large-scale  network,  thus  a  distributed 
cooperative  routing  and  resource  allocation  algorithm  is  needed  to  take  link  achievable  rate,  link 
interference  to  the  PU  receiver,  and  mutual-information  accumulation  into  consideration.  In  fact, 
distributed  algorithms  have  been  studied  in  [Draper20 11  Cooperative]  and 
[Urgaonkar2012Optimal],  but  they  do  not  take  the  interference  from  SUs  to  PUs.  Although 
[Chowdury2011CRP],  [Caleffi2012OPERA]  and  [Chen2012A]  consider  the  interference  power 
constraint  of  PUs,  they  cannot  be  directly  used  in  CRNs  with  mutual-information  accumulation. 
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Figure  5:  A  simple  three-node  SU  network. 

Take  a  simple  three-node  SU  network  as  an  example.  As  shown  in  Fig.  5,  assume  that  node  1 
and  node  2  have  already  successfully  decoded  the  message  while  node  3  still  needs  to  receive  Bq 
bits  of  mutual  information  before  successfully  decoding  the  message.  C12  ,  C13  ,  C23  are 
respectively  the  link  achievable  rates  among  node  1,  2  and  3. 1 j  is  the  interference  threshold  set 
by  the  PU  receiver  and  the  total  bandwidth  of  the  subband  is  normalized  to  Wj  =  1.  and  62  are 
the  spectrum  bandwidth  allocated  to  node  1  and  node  2.  The  interference  power  spectral  density 
of  node  1  and  node  2  to  the  PU  receiver  are  thus  respectively  =  h^rPs  ^  and  1 2  =  h2rPs^  then  the 
optimal  spectrum  resource  allocation  can  be  expressed  as 
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(31) 


The  01  and  0|  can  be  applied  to  minimize  the  end-to-end  delay.  From  the  equation,  we  find 
the  delay-optimal  cooperative  route  for  underlay  CRNs  depends  on  both  the  interference  condition 
and  the  link  achievable  rate  condition;  and  under  some  interference  and  link  achievable  rate 
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conditions,  both  node  1  and  node  2  have  the  opportunity  to  participate  in  the  transmission 
simultaneously  (see  the  fifth  and  sixth  items). 

The  above  equation  is  the  optimal  solution  for  a  three-node  network.  To  extend  the  result  to 
multiple  nodes  scenario  the  problem  can  be  solved  by  introducing  the  concept  of  a  virtual  node.  If 
there  are  two  nodes,  say  node  1  and  node  2,  doing  concurrent  transmission  at  the  same  slot,  the 
two  nodes  combined  can  be  regarded  as  a  virtual  node  2.  If  a  node  decodes  the  packet  and 
determines  that  a  virtual  node  is  transmitting  in  the  network,  this  node  will  not  transmit 
concurrently  with  the  virtual  node,  and  within  a  slot,  there  are  at  most  two  SUs  sharing  PU’s 
spectrum  and  transmitting  concurrently.  Further,  for  the  underlay  CRN  with  M  primary  pairs, 
within  a  slot  there  are  at  most  M  +  1  SUs  sharing  PU’s  spectrum  and  transmitting  concurrently. 

The  following  algorithm  determines  the  distributed  routing  and  resource  allocation  algorithm 
for  underlay  CRNs  with  mutual-information  accumulation.  The  computational  complexity  of  this 
distributed  algorithm  for  each  node  is  0(1),  and  the  overhead  of  our  distributed  algorithm  can  be 
quite  small. 


Algorithm  2:  Distributed  routing  and  resource  allocation 
algorithm  for  underlay  CRNs  with  MIA 

1  Initialization:  (For  node  i)  Estimate  its  interference 
power  spectral  density  li  to  the  PU  receiver;  estimate  its 
channel  to  the  destination  and  calculate  CiNl  obtain  the 
information  of  current  transmitting  node  li-i  and 
C(i-i)N  from  the  header  of  the  data  packet  it  has 
decoded; 

2  calculate  0i,02  by  equation  (7); 

3  if  01  =  0,  02  >  0  then 

4  Set  01  =  0  in  the  control  packet  and  send  it  out  to 
inform  current  transmitting  node  to  stop  transmitting; 
set  V_flag  to  be  FALSE;  add  CiN,  U  in  the  header 
and  start  to  transmit; 

5  else  if  01  >  0,  02  >  0  then 

6  Set  01  in  the  control  packet  using  the  computed  value 
of  01  in  step  2;  send  out  the  control  packet;  in  the 
header  of  a  data  packet,  set  V_flag  to  be  TRUE;  set 
CiN  to  be  the  equivalent  achievable  rate  of  the 
virtual  node;  set  li  to  be  It',  transmit  data  packet 
with  bandwidth  02; 

7  else 

8  node  i  keeps  silent. 
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3.4  USRP  Demo 


USRP  N210s  are  used  to  demonstrate  our  results.  More  specifically,  we  construct  a  five-node 
DSA/CRN  network  and  using  GNU  radio  to  implement  the  protocol.  Two  applications  are 
demonstrated:  1)  image  file  transfer  over  the  SU  network,  and  2)  video  file  transfer  over  the  SU 
network. 

The  GNU  radio  node  design  can  be  found  in  the  following  figure. 


The  GNU  physical  layer  design  can  be  found  in  the  following  figure. 
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Figure  7:  PHY  Layer  Design  (GNU  Radio) 


RESULTS  AND  DISCUSSION 

Spectrum  Sensing  for  Dynamic  Spectrum  Access  Networks 


Figure  8:  Optimal  sensing  time  versus  total  transmit  power  of  SU. 
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Fig.  8  shows  the  relationship  between  SU’s  optimal  sensing  time  and  SU’s  total  transmit  power. 
Eaeh  point  in  the  figure  is  the  average  results  of  5000  simulation  runs  and  P(//o)  =  0.7.  We  can 
see  that  when  PU  has  stringent  rate  loss  constraint,  as  SU’s  power  budget  increases,  SU  will  spend 
more  time  on  sensing.  However,  if  PU  can  tolerant  much  larger  rate  loss,  as  SU’s  power  budget 
increases,  SU’s  optimal  sensing  time  decreases.  This  difference  is  caused  by  the  exchange  of  the 
dominant  constraint  of  the  system.  Another  interesting  finding  is  that  when  SU’s  power  budget  is 
low,  the  optimal  sensing  time  of  SBSS  will  be  smaller  than  that  of  spectrum  overlay.  But  as  SU’s 
power  budget  increases,  SU’s  optimal  sensing  time  in  SBSS  will  be  larger  than  that  of  spectrum 
overlay. 


Figure  9:  SU’s  throughput  versus  SU’s  sensing  time 

Fig.  9  compares  the  performance  of  SBSS  and  spectrum  overlay  in  multi-carrier  case  under 
different  sensing  times.  The  black  line  represent  the  throughput  of  cognitive  system  with  only 
power  constraint.  The  total  number  of  carriers  is  set  as  30,  PU’s  power  is  25  dBm  and  SU’s  total 
power  budget  is  20  dBm  (these  two  parameters  are  set  as  30  times  of  corresponding  values  in  the 
single-carrier  case).  PU’s  rate  loss  6  is  set  as  0.1.  From  this  picture,  we  can  see  that  in  multi-carrier 
case,  SBSS  can  improve  the  spectrum  efficiency  one  step  further  compared  to  the  corresponding 
single-carrier  case,  specifically,  under  the  optimal  sensing  time  (0.2ms),  SBSS  will  achieve  15.2% 
more  spectrum  efficiency  than  spectrum  overlay. 

4.2  Rateless  Coding  based  Cooperative  Dynamic  Spectrum  Access  Networks 
4.2.1  Improving  the  throughput  of  CRN  with  mutual-information  accumulation 
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A 

Figure  10:  Maximum  stable  throughput  p_S  versus  the  throughput  selected  by  primary  node  X_P. 


Fig.  10  shows  the  the  maximum  stable  throughput  Hs  versus  the  throughput  selected  by  the 
primary  node  Ap.  We  can  see  from  the  figure  that  over  any  Ap  <  exp  =  0.37,  EA  and  mutual- 

information  accumulation  achieve  higher  stable  throughput  of  the  cognitive  link  than  that  without 
any  accumulation,  with  mutual-information  accumulation  protocol  higher  than  EA  protocol.  That 
is  due  to  increasing  the  probability  of  successful  primary  packets  reception  at  PR,  thereby  creating 
more  transmission  opportunities  for  ST. 


Figure  11:  Maximum  throughput  of  the  secondary  user  fis  versus  Ysp  (Ap  «  exp 


(&).0.37,. 


Fig.  11  plots  the  maximum  stable  throughput  Hs  of  the  secondary  link  for  varying  values  of 
Ysp-  As  Ysp  increases,  the  relaying  gain  increases  for  all  the  scenarios  under  study,  since  most  of 
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the  traffic  is  redirected  to  the  secondary  link  due  to  the  better  channel  conditions  on  ST-PR  link. 
Moreover,  we  can  observe  that  the  maximum  stable  throughput  Hs  for  EA  and  mutual-information 
accumulation  protocols  achieve  higher  throughput  than  that  without  any  accumulation.  In  addition, 
the  performance  of  mutual-information  accumulation  is  slightly  higher  than  that  of  EA  when  Ysp 
is  low,  and  the  two  become  similar  in  performance  as  Ysp  becomes  high. 

4.2.2  Cooperative  Routing  for  Underlay  CRNs  using  Mutual-Information  Accumulation 

Fig.  12  gives  an  example  to  compare  the  routing  and  resource  allocation  solutions  between 
traditional  multi-hop  delay-optimal  routing  without  mutual-information  accumulation  and 
cooperative  routing  using  mutual-information  accumulation  in  underlay  CRNs.  All  the  nodes  are 
distributed  in  the  area  of  1  X  1  (meter  xmeter).  The  PUT  is  placed  at  [0.4, 0.8]  and  the  PUR  is 
located  at  [0.6, 0.8].  There  are  altogether  16  SU  nodes  in  the  network:  the  SU  source  (node  1)  is 
placed  at  [0.1, 0.4]  and  the  SU  destination  (node  16)  is  located  at  [0.9, 0.4].  The  other  14  SU  relay 
nodes  are  uniformly  distributed  in  the  area.  The  total  system  bandwidth  is  set  to  be  Wj  =  1  Hz 
and  the  total  system  energy  is  set  to  be  Ej  =  20  Joule.  As  shown  in  The  corresponding  bandwidth 
allocation  is  labeled  in  the  figure.  For  example,  when  the  source  node  transmits  to  node  15,  it  will 
use  82%  of  the  total  bandwidth  due  to  the  interference  power  constraint.  The  route  for  cooperative 
routing  with  mutual-information  accumulation  is  [1,  (8,1),  (8,3),  (8,6),  15, 4,  (4,2),  16],  where 
the  item  of  (8, 1)  means  that  node  8  and  node  1  can  be  regarded  as  a  virtual  node  to  share  PU’s 
spectrum  jointly  in  the  second  slot.  Within  this  virtual  node,  node  1  will  occupy  54.3%  of  the  PU’s 
spectrum  while  node  8  will  occupy  the  rest  45.6%  of  the  spectrum.  From  Fig.  12,  we  can  see  that 
the  delay-optimal  route  of  cooperative  routing  using  mutual-information  accumulation  for 
underlay  CRN  is  very  different  from  that  using  traditional  multi-hop  routing.  It  is  important  to  note 
that  the  cooperative  routing  and  resource  allocation  solution  in  [Draper20 11  Cooperative]  can  be 
regarded  as  a  special  case  of  our  algorithm  by  setting  Ij  ^  co. 
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Figure  12:  Sample  underlay  CRN.  Delay-optimal  cooperative  routing  using  Mutual-Information  Accumulation  (solid, 
Mutual-Information  Accumulation  Routing)  and  traditional  multi-hop  delay-optimal  routing  (dashed,  TM  Routing). 


4.3  Distributed  Resource  Allocation 

4.3.1  Cooperative  Retransmission  for  massive  MTC 


Figure  13:  Outage  probability  as  a  function  of  cooperative  region. 

Fig.  13  shows  simulation  results  as  well  as  numerical  results  on  the  probability  of  outage  for 
both  conventional  retransmission  P’"®  and  the  proposed  cooperative  strategy  P^°.  In  this  figure, 
each  point  is  an  average  over  5000  realizations  and  A  =  10“^m“^,  p  =  0.005,  a  =  4,  and  z  = 
400  m.  From  Fig.  13,  one  can  see  that  the  introduced  cooperative  strategy  can  significantly  reduce 
the  outage  probability  compared  to  conventional  retransmission.  Specifically,  with  a  decoding 
threshold  of  2  dB,  the  outage  probability  of  massive  MTC  network  with  cooperative  region  of  65 
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m  is  14.5%  while  corresponding  conventional  retransmission  has  an  outage  probability  of  39.7%. 
Note  that  the  outage  probability  of  conventional  retransmission  is  higher  than  the  outage 
probability  of  the  first  time-slot  transmission  (which  is  37.0%)  due  to  the  temporally  correlated 
interference.  The  reason  why  our  introduced  cooperative  strategy  can  reduce  the  outage  probability 
is  because:  1)  the  relay  is  selected  only  when  it  can  decode  the  packet;  2)  the  selected  relay  within 
the  cooperative  region  is  much  closer  to  the  BS.  Furthermore,  we  can  find  that  as  the  radius  of 
cooperative  region  reduces,  the  outage  probability  of  cooperative  communication  will  approach  to 
the  outage  probability  of  conventional  retransmission.  This  is  because  as  the  cooperative  region 
reduces,  the  probability  of  finding  potential  relays  reduces. 


4.3.2  Distribute  Cooperative  Routing  and  Resource  Allocation  in  Underlay  CRNs 


Figure  14:  Delay  distribution:  traditional  multi-hop  delay-optimal  routing,  distributed  and  centralized  solution,  optimal 

solution  obtained  through  exhaustive  search. 

Fig.  14  shows  the  cumulative  delay  distribution  for  the  following  four  approaches:  traditional 
multi-hop  delay  optimal  routing,  centralized  algorithm  for  mutual-information  accumulation, 
distributed  algorithm  for  mutual-information  accumulation,  and  optimal  algorithm  for  mutual- 
information  accumulation  obtained  by  exhaustive  search.  6  SU  nodes  are  randomly  deployed  in 
the  network  as  SU  relays.  Ij  =  3W  and  a  =  3.  The  average  delay  of  the  traditional  multi -hop 
delay-optimal  routing  is  11.546  sec,  while  the  average  delay  of  centralized  and  distributed 
cooperative  routing  using  mutual-information  accumulation  are  1.620  sec  and  1.622  sec, 
respectively.  The  optimal  delay  obtained  by  exhaustive  search  is  1.608  sec.  On  average,  both  the 
centralized  and  distributed  cooperative  routing  algorithms  can  reduce  up  to  85.9%  of  the  end-to- 
end  delay  compared  to  traditional  multi-hop  delay-optimal  routing.  Furthermore,  compared  to  the 
optimal  route  and  resource  allocation  strategies  using  exhaustive  search,  our  proposed  centralized 
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and  distributed  algorithms  only  add  0.75%  and  0.87%  additional  delay,  respectively.  Under  the 
same  simulation  parameters,  if  we  change  the  number  of  SU  nodes  to  48,  the  centralized  and 
distributed  algorithms  can,  on  average,  reduce  77.7%  and  77.5%  of  the  end-to-end  delay  compared 
to  the  traditional  multi-hop  delay-optimal  routing.  Finally,  this  simulation  figure  indicates  that  the 
introduced  distributed  algorithm  can  obtain  almost  the  same  performance  on  end-to-end  delay 
compared  to  the  centralized  algorithm. 


4.4  USRP  Demo  and  Graphical  User  Interface 

In  June  2015,  cooperative  routing  using  rateless  coding  for  a  three-node  network  was 
demonstrated  at  the  Griffiss  Institute  using  USRP  N210s. 

Demo  Setup:  The  details  of  the  demo  are  as  follows:  1)  2  laptops  controlling  3  USRP  N210s;  2) 
3  USRP  N210s  serving  as  Source,  Relay,  and  Destination;  3)  carrier  frequency  is  2.4  GHz;  4) 
network  bandwidth  is  1.25  MHz;  5)  sampling  rate  is  5  Mbps;  6)  slot  duration  for  TDMA  is  100 
ms;  and  7)  size  of  Raptor  symbol  is  2056  bytes  with  Raptor  K  =  25.  This  setup  can  be  seen  most 
clearly  in  the  following  figure. 


Laptop  1 

Figure  15:  Network  Setup  for  USRP  Demo 


Medium  Access  Control  (MAC)  layer:  The  frame  structure  can  be  seen  most  clearly  in  the 
following  figure.  It  is  important  to  note  that  the  MAC  layer  control  signaling  is  implemented  over 
the  air  instead  of  relying  the  wired  line  communication  to  serve  as  the  control.  Due  to  the  wireless 
control,  we  need  to  reserve  enough  system  resource  for  wireless  ACK/NAK  to  reach  the  receiver 
from  the  transmitter. 
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Figure  16:  MAC  Frame  Structure 


Graphical  User  Interface  (GUI):  GUI  design  was  implemented  and  demonstrated  in  May  2016 
during  the  TIM  meeting  with  AFRL.  The  GUI  utilizes  remote  database  (MYSQL)  and  shows  real¬ 
time  accumulated  mutual-information  from  each  transmission.  This  gives  a  visual  representation 
of  the  benefits  of  rateless  coding-based  cooperative  routing.  The  developed  GUI  can  be  seen  more 
clearly  in  the  following  figure. 


Figure  17:  GUI  for  Cooperative  Routing 


Demo  Outcome:  In  both  demonstrations  (June  2015  and  May  2016),  we  showed  the  following 
features  in  different  communication  layers:  1)  Modulation  and  demodulation,  as  well  as  rateless 
encoder/decoder  (RFC5053)  at  the  physical  layer;  2)  Time  division  multiple  access  (TDMA)  and 
routing  in  the  MAC  and  Network  layer;  3)  File/image  file  transfer,  live  video  streaming  and  GUF 
remote  database  control  in  the  application  layer. 

During  the  USRP  demo  we  showed  that  the  transmission  delay/time  for  non-cooperative 
communication  (from  the  source  to  the  destination  directly)  is  57  seconds  while  the  transmission 
delay/time  for  the  cooperative  strategy  we  developed  is  30  seconds.  This  results  a  cooperative  gain 
of  1.9  which  matches  very  well  with  our  analytical  characterization  (2.0). 
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5.0  CONCLUSIONS 

5.1  Spectrum  Sensing  for  Dynamic  Spectrum  Access  Networks 

We  studied  optimal  sensing  time  and  power  allocation  strategies  for  sensing  based  spectrum 
sharing  networks  subject  to  a  PU’s  rate  loss  constraint  (RLC)  and  power  constraint.  Specifically, 
in  the  single-carrier  case,  we  compared  the  performance  (SU’s  optimal  throughput  and  optimal 
sensing  time)  of  SBSS  with  that  of  spectrum  overlay.  In  the  multi-carrier  case,  we  proposed  an 
efficient  algorithm  to  find  the  optimal  sensing  time  and  power  allocation.  In  single-carrier  case, 
SBSS  achieved  more  throughput  than  spectrum  overlay  when  SUT  is  far  away  from  PUR  or  PU 
can  tolerant  large  rate  loss;  in  multi-carrier  case,  SBSS  achieved  higher  spectrum  efficiency  than 
the  corresponding  single-carrier  case  by  exploiting  additional  power  allocation  flexibility. 

5.2  Rateless  Coding  based  Cooperative  Dynamic  Spectrum  Access  Networks 

Mutual-information  accumulation  is  an  enhanced  technique  to  increase  the  performance  of  the 
dynamic  spectrum  access  network.  The  receiver  can  decode  its  packet  data  if  the  mutual- 
information  accumulation  exceeds  a  threshold  via  multiple  transmissions  from  multiple 
transmitters,  and  every  transmission  is  useful  to  contribute  some  mutual  information.  In  practice, 
mutual-information  accumulation  can  be  realized  through  rateless  codes  in  which  the  transmitter 
generates  unlimited  number  of  encoded  packets,  such  that  the  receiver  can  decode  the  data  after 
receiving  a  sufficient  number  of  encoded  packets,  irrespective  of  which  ones  it  has  received.  This 
improves  the  likelihood  and  ability  of  multiple  nodes  to  cooperatively  serve  users  and  improve 
DSA  network  performance. 

From  the  description  of  the  mutual-information  accumulation  used  in  CRN  to  improve  the 
throughput  of  the  network,  and  minimize  the  end-to-end  delay  in  multi-hop  underlay  CRN,  we 
found  that  mutual-information  accumulation  outperforms  the  others  significantly.  In  specific, 
applying  mutual-information  accumulation  in  CRN,  where  SUs  can  decode  and  relay  the  packets 
of  PUs,  performs  best  compared  with  energy  accumulation  and  without  accumulation  scenarios. 
Also,  mutual-information  accumulation  was  used  to  find  optimal  routing  and  resource  allocation 
scheme  to  minimize  the  end-to-end  delay  of  the  second  users  while  the  QoS  of  primary  users  are 
guaranteed.  The  end-to-end  delay  was  reduced  up  to  85.9%  compared  with  the  traditional  scheme 
without  mutual-information  accumulation  [Chen20 15 Cooperative]. 


5.3  Distribute  Resource  Allocation 

Centrally  determining  resource  allocation  for  all  candidate  transmission  nodes  in  a  large-scale 
network,  if  even  possible,  is  an  exceedingly  expensive  operation.  A  potential  solution  is  to  perform 
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distributed  resource  allocation  in  which  the  resource  allocation  is  determined  locally.  In  a 
distributed  resource  allocation  scheme,  a  node  does  not  need  to  know  the  information  of  the  whole 
network,  instead  it  makes  a  decision  based  only  on  the  information  of  its  own  and  surrounding 
nodes.  Also,  if  one  node  is  out  of  charge,  then  a  new  node  can  be  found  locally,  and  quickly. 
Compared  with  the  centralized  resource  allocation,  the  distributed  scheme  is  a  quite  cost-efficient, 
flexible,  and  robust,  and  is  widely  used  in  many  research  fields  such  as  massive  MTC  and 
cooperative  routing  in  multi-hop  CRNs. 

In  MTC  networks,  the  inactive  MTDs  were  allowed  to  act  as  relays  for  outage  MTDs.  Using 
point  process  theory,  appropriate  relay  MTDs  were  selected  to  minimize  the  outage  probability  of 
this  kind  of  massive  MTC  network.  Both  numerical  and  simulation  results  demonstrated  the  great 
promise  of  cooperation  of  the  relay  MTD  and  the  source  MTD  in  massive  MTC  networks. 

In  cooperative  routing  multi-hop  CRNs,  an  efficient  distributed  routing  and  resource 
allocation  scheme  was  proposed,  and  it  was  shown  that  the  complexity  of  the  distributed  scheme 
was  lower  than  the  centralized  scheme.  Simulation  results  showed  that  the  performance  of  the 
distributed  scheme  is  very  close  to  the  centralized  scheme. 

5.4  USRP  Hardware  Demo 

Hardware  demonstrations  using  USRP  N210s  were  conducted  during  two  TIMs  in  AFRL.  The 
basic  features  of  the  rateless  coding-based  cooperative  routing  were  implemented  in  GNU  radio 
and  the  protocols  in  almost  all  communication  layers  (from  PHY  to  APP)  were  changed.  The  demo 
showed  that  performance  gains  of  rateless  coding-based  cooperative  routing  were  realistic  and 
achieved  in  practice.  Furthermore,  the  achieved  performance  gain  matched  very  well  to  theoretical 
results. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


ARQ 

automatic  repeat  request 

BS 

base  station 

CRN 

cognitive  radio  network 

CSI 

channel  state  information 

DAS 

dynamic  spectrum  access 

D2D 

device-to-device 

EA 

energy  accumulation 

LP 

linear  programming 

MTC 

machine  type  communication 

MTD 

machine  type  device 

PST 

power  spectrum  density 

PT 

primary  transmitter 

PR 

primary  receiver 

PU 

primary  user 

PUR 

primary  user  receiver 

PUT 

primary  user  transmitter 

QoS 

quality  of  service 

RLC 

rate  loss  constraint 

SBSS 

sensing  based  spectrum  sharing 

ST 

secondary  transmitter 

SR 

secondary  receiver 

SU 

secondary  user 
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SUR 

SUT 


secondary  user  receiver 
secondary  user  transmitter 
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